Efficient error characterization in Quantum Information Processing 
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We describe how to use the fidelity decay as a tool to characterize the errors affecting a quantum 
information processor through a noise generator G T . For weak noise, the initial decay rate of the 
fidelity proves to be a simple way to measure the magnitude of the different terms in G T . When the 
generator has only terms associated with few-body couplings, our proposal is scalable. We present 
the explicit protocol for estimating the magnitude of the noise generators when the noise consists of 
only one and two-body terms, and describe a method for measuring the parameters of more general 
noise models. The protocol focuses on obtaining the magnitude with which these terms affect the 
system during a time step of length r; measurement of this information has critical implications for 
assesing the scalability of fault-tolerant quantum computation in any physical setup. 



I. INTRODUCTION 

One of the biggest challenges in the physical realization 
of quantum information processing (QIP) is the precise 
control of the system. In order to achieve this, we must 
be able to characterize the noise sources affecting exper- 
imental setups. A first answer to this problem was given 
by quantum state and quantum process tomographics 
(QST, QPT), which allow one to reconstruct the dynam- 
ics occurring in the probed system during a given time 
step But beyond a few qubits, a complete characteri- 
zation of the errors through this method is not practical, 
since for a system of n qubits it would require 0(2 2n ) 
measurements |2j. 

This has motivated the search for protocols showing 
particular features of the system dynamics but with a 
scalable implementation, as in 

S II lH Recently, 
schemes using random maps have received special atten- 
tion: just as random numbers play a fundamental role in 
classical information theory, random unitary operators 
and random quantum states are also very useful com- 
ponents for quantum information theory and processing 
Q. Exact random operators are exponentially hard to 
implement, turning any of these proposals into unscal- 
able ones. However, it was recently suggested |8| that 
pseudo-random operators displaying the principal fea- 
tures of random behavior can be constructed efficiently. 

Following this, efforts have been made to find sig- 
natures of the environment in evolutions generated by 
random operators. Emerson et al. [t| have obtained 
results linking the fidelity decay with the strength of the 
noise affecting the system. 

Here we describe how to extract additional informa- 
tion about the errors. In our approach the noise over 
a time interval t is represented by an error operator 
E = cxp(— iG T ), where G T is the noise generator. The 
composition of G T gives rise to a variety of noise models. 
These models with their assumptions are stated in Sec. 
II and further illustrated in Appendix A. Several results 
were obtained showing an exponential-like fidelity decay, 



which we report at the end of Sec. II and further de- 
velop in Appendix B. Moreover, we derived an analytical 
expression for the initial decay rate, under the assump- 
tions of weak noise and relatively good control. As we 
will show, this particular result allows us to determine 
the magnitude of the terms in G T with different Ham- 
ming weights for a given subset of the n qubits. (The 
Hamming weight of a product operator -i.e., an opera- 
tor that is a product of the Pauli matrices a x , a y , u z 
or the Identity operator I for each qubit- is the number 
of factors which are different from I.) These results are 
presented in Sec. Ill, with some more general formulae 
deferred to Appendix C. Sec. IV is devoted to practical 
implementation issues. Finally, in Sec. V, we present our 
conclusions and an outline of future work. 



II. DESCRIPTION OF THE NOISE MODEL 
AND GENERAL RESULTS 

We start with the standard definition of fidelity after 
motion reversal for a pure state (also called Loschmidt 
Echo), 



/(<) = |(*„ | U\t)U P (t) | * > | 



(1) 



where U is the "perfect" evolution of the system, while 
Up represents a perturbation of U. As it has been shown 
before (see [j, [l(| and references therein) , the fidelity af- 
ter motion reversal encodes information about the perfect 
evolution and/or about the perturbation. Our aim is to 
obtain information about the perturbation, which for us 
consists in the errors resulting from uncorrelated imper- 
fect implementation of gates and environmental/external 
effects. The question is how to choose the perfect evolu- 
tion in order to only reveal information about the errors. 
Random operators are a reasonable choice, since random- 
ness enables the exploration of the whole Hilbcrt space 
while removing the effect of irrelevant parameters (such 
as particular choices of initial states). In addition, the 
use of random operators suppress the effect of possible 
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FIG. 1: Circuit representation of the algorithm. The opera- 
tors R are randomly drawn in every step. The errors E are 
also allowed to vary in each step, (a) The algorithm seen as a 
fidelity decay scheme, as given by eqs. @ & 0. (b) Equiva- 
lent algorithm, seen in a twirling fashion. The equivalence of 
both designs is due to the average over the random operators 
R which are taken to be invariant under the Haar measure. 



time-correlations between the terms in the noise genera- 
tor. 

To incorporate initial states other than pure states, we 
took a generalization of for density matrices 



/(*) 

with p mr (t) 



Tr[p mr (t)po] 
R\ . . . R\ E t Rt 
x R\ E\ . 
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where po is the initial state of the system, the operator 
E t represents the errors cumulated during one step of the 
algorithm at time t and R t represents some random op- 
erator in the Hilbert space Jtiffo of n qubits (N — 2"). 
The subscript mr stands for motion reversed. The algo- 
rithm is depicted in Fig. Q (a). We take the operators 
Et to be unitary. Since we use random operators, we 
are actually concerned with the average fidelity (fit))' 
the brackets denote averaging over the considered group 
of random operators, which we take to be invariant un- 
der the Haar measure. This average makes this fidelity 
decay scheme equivalent to a twirling scheme [l 11 , as de- 
picted in Fig. ^ (b). This equivalence can be derived 
from the unitary invariance of the Haar measure. Being 
a stcp-by-step twirling of E, our proposal breaks up the 
action of E making the state at one step depend only on 
the previous one. Therefore, we can already conjecture 
an exponential-like decay for (fit)), as in the work by 
Emerson et al. Q. 

Defining the fidelity between two mixed states as in 
PJl may seem arbitrary; there is no unique generalization 
of for initial states that are not pure (see for example 
It will become clear later that this mathematical 
entity reflects the quantity we measure in order to 
implement our proposal, and for convenience we shall 
call it fidelity throughout this work. 

We chose the random operators to be R s = R^ ® 
Ri^ <g> . . . Rs , where Rg is a random rotation of the 
qubit j (the resulting algorithm is depicted on Fig. |5J). 
Another possible choice, consisting of R s being uniformly 
drawn from U(N), has already been studied in this 
led to general and closed results for (f(t)), essentially 
showing a universal exponential decay depending only 
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FIG. 2: Circuit representation of the algorithm, eqs. (2) & 
©, choosing individual rotations as the random operators. 



on the purity of the initial state (that is Tr[pg]) and on 
the global strength of the noise, quantified by the trace of 
the superoperator describing the non-unitary dynamics. 
Although this is a useful analytical result, this strong 
randomization scheme doesn't yield any information on 
the noise structure, hinting at the usefulness of a weaker 
form of randomization. 

The Rg^ are drawn uniformly from SU (2) with respect 



to the invariant Haar measure, 
computational basis is [l3| 



Their expression in the 



RP = 



,)e lXl ' s cos(4>j, s )e 



with ip and \ drawn uniformly from the interval [0, 2tt), 
and cj) = arcsin (V?) with £ uniformly distributed in 
[0,1). 

Notice that these random operators can be efficiently 
implemented, since they are single-qubit operations 
with a suitable gate decomposition fl4l |. In this respect, 
we're not affected by the efficiency issue that arises 
in the implementation of random operators in U(N) 
with arbitrarily large N, which leads to the use of 
pseudo-random operators. 

We assume we have relatively good control of the sys- 
tem under study, so the random rotations can be imple- 
mented with sufficient accuracy and then the errors are 
only present in the error operators Et- This is a reason- 
able hypothesis since we can always make the magnitude 
of the errors stemming from E relatively larger than the 
ones in the random rotations by increasing the implemen- 
tation time of E. A fair exploration of the errors affecting 
the system can be achieved by trying to implement the 
identity operator I, i.e., trying to prevent any evolution 
of the system for a certain time r. 

Thus we consider £ to be a deviation from I in the 
form 



JV -l 

E = exp(-iG T ) = cxp ( -i ^ xiO, 

i=i 



(4) 



where G T is the generator of E. We can regard E to be 
a residual operator resulting from the action of the noise 
during a time r through an effective Hamiltonian HG t /t, 
for the time step under consideration. Without losing 
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generality, the generator is decomposed in the product 
operator basis Of. 

n 

0/ = (g)Op ) (5) 

j=i 

where each 0\ is an operator in the space of qubit j 
and it is either a Pauli matrix or the Identity (thus the 
Oi are hermitian), but at least one factor in each Oi is 
a Pauli matrix (thus the Oi are traceless). For the error 
operators E to be unitary, the coefficients xi must be 
real numbers. 

We can distinguish three classes of effective noise mod- 
els depending on the time variation of E: 

C - Coherent: E remains the same at all times. 

IL - Incoherent with long correlation time: E remains 
the same during approximately the time required 
to implement one realization of f(t). 

IS - Incoherent with short correlation time: E changes 
from gate to gate. 

In types IL and IS, the change of E will be given by a 
change of the coefficients xi i n G T ; these are randomly 
drawn according to a given distribution f ({%/})■ Co- 
herent noise gives unitary errors. Incoherent scenarios 
give rise to non-unitary errors; we call them "incoher- 
ent" since the superoperator arising in these cases has 
the form of an (incoherent) average over the parameters 
characterizing the unitary operation: 

S B (p) = J E ({xi}) pE^ixi}) P({xi}) d X i (6) 

Notice that in each realization, although at each step 
the action of the error operator is unitary, its variation 
in time introduces a net non-unitary operation when 
averaging different realizations to obtain (f(t)). A brief 
illustration of these models is given in Appendix A for 
the one qubit case. 

We expect that some of all the terms we allow in the 
noise generator will be negligible. We will analyze G T 
for truncated sums over multi-body terms Oi up to a 
given Hamming weight. For this we introduce a more 
specific labeling of the terms in G T . When necessary, we 
shall denote the xi as Xj'k > where j, k, . . . label qubits, 
and together with p, q, . . . = {x, y, z}, they indicate 
that that particular term is a product of the Pauli 
matrices a p for qubit j, a q for qubit k, etc., and the rest 
of the factors are just the Identity for the other qubits. 
Therefore, the one-body terms (Hamming weight 1) go 
with coefficients Xjt Xji Xji two-body terms (Hamming 
weight 2) are X P \k> e ^ c - To avoid double counting of 
multi-body terms, the labeling of the qubits must obey 
j < k < . . . and so on. 



As we already mentioned, the non-negligible coeffi- 
cients can in general be drawn from any given distribu- 
tion P({xi})- We studied (/(£)) for some specific cases, 
including only one-body terms (analytical and numerical 
results) and one-body and two-body terms (numerically) . 
The two distributions considered were: 1) a constant dis- 
tribution (x P j = oi, x P % = P V j, p, k, q); 2) each x P 
(resp. Xjk) randomly drawn from a Gaussian distribu- 
tion with mean value a ((3) and standard deviation a a 
(erg)- We will refer to the coefficients xi or to a, /3, a a , 
ap collectively as the "noise strength" x, an d the pow- 
ers of x wm include any monomial combination of degree 
equal to the given power. 

Numerical calculations of (/(t)) show an exponential- 
like decay. In particular we observe: 

i) Linear initial decay: (/(t)) ~ /o(l — 7*) for t suffi- 
ciently small and with /o = Tr^] (we shall call 7 
the initial decay rate); 

ii) Constant long-time limit: (/(*)) — > 1/N for t — > 00. 

The scales u t sufficiently small" or t — ► 00 are set by 
the strength of the noise x- O ur numerical calculations 
ranged up to a strength of 0.4. For higher strength, the 
saturation value 1/N is reached in only a few steps and 
not much can be extracted from this fidelity decay. An 
example to illustrate these calculations is given in Fig. 
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FIG. 3: Examples of the numerical calculation of (f(i)} with 
100 realizations for 8 qubits, initially all in the |0) state. Here 
we took G T with one-body terms only (□, A. O) and with one- 
body and all the two-body terms (*, •, 0). 
□ : type IL, P Gaussian with a = f3 = 0, a a = 0.05, up — 0. 

A type C, P constant with a = 0.05, ,3 = 0. 

O: type IS, P Gaussian with a = f3 = 0, a a — 0.05, Gp = 0. 

*: type IL, P Gaussian with a = j3 = 0, a a — erg = 0.05. 

• : type C, P constant with a = j3 = 0.05. 

0: type IS, P Gaussian with a = j3 = 0, a a — 073 = 0.05. 



Analytical expressions of (f(t)) can be obtained if we 
consider one-body terms only and a separable initial 
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state. Our closed results are in exact agreement with 
the behavior described above. The slope of the decay is 
of order 0(x 2 ) and is the same for the different types of 
noise C, IL and IS. The analytical results were derived 
with mathematical tools developed in ^(|, a good pre- 
sentation thereof can be found in [T^ . This approach has 
been used in 0] to study the fidelity decay of perturbed 
quantum chaotic maps. 

We defer further details about the calculations for the 
full (f(t)), since the most interesting results arise from 
our studies of the initial decay rate of (f(t)). The reader 
is referred to Appendix B for a more complete report on 
the former. 

To conclude this Section we revisit the following point: 
our work indicated that the fidelity decay after motion 
reversal is initially linear in time, also in agreement 
with the results already published by Emerson et al. . 
While this seems to contradict previous well established 
results reporting a universal quadratic decay [Tof. these 
two statements do not in fact contradict each other, since 
the random dynamics studied here is not considered by 
previous work, which presumes the use of a constant 
evolution operator. In our case, the evolution is given 
by random rotations which vary in each step and in 
each realization, and we then study the evolution of 
an ensemble-averaged state. Even when some relations 
between our error operator E and the perturbation 
can be drawn (which originally motivated the use of a 
fidelity decay scheme), the nature of the calculation is 
different. 



III. THE INITIAL DECAY RATE 7 

Numerical and analytical evidence supports the con- 
jecture that the initial decay is linear in t, with an initial 
decay rate 7. Moreover, this initial decay rate depends 
only on the noise strength (the magnitude of the xi) an d 
not on the particular time variation of E. This can be 
seen in Fig. where the initial decay is the same for the 
different types C, IL and IS (described in Sec. II) as long 
as the general noise strength is the same. In addition, 
the dependence on the strength is quadratic (Fig. |3J). 

The first clear evidence comes from the initial slope of 
the analytical expressions we obtained; this result is of 
course limited to noise with only one-body terms. For 
noise including two-body terms, we fitted the initial de- 
cays (as shown in Fig. [3J, obtaining a quadratic depen- 
dence on the governing parameter. 

We see then that the initial decay rate is a measure of 
the noise strength independently of the type of noise. By 
defining precisely the initial decay rate 7 in terms of the 
fidelity after the first step, 




Noise strength % 

FIG. 4: Values of the initial decay rate 7 obtained with a 
linear fitting of (/(£)) in the small t regime. In practice, this 
linear regime was given by all the points t with {/(£)) > fu m . 
Here fu m = 0.9 and n = 8. Solid lines: the quadratic fitting. 
O: type C, P constant with a — j3 = \\ au the terms in G T . 
• : type IL, P Gaussian with a — /3 = 0, a a — up — \\ au the 
terms in G T . 

*: idem O, but only one-body and first-neighbor terms. 
□ : idem •, but only one-body and first-neighbor terms. 
<0>: idem O, but only one-body terms. 
A: idem •, but only one-body terms. 



we can actually obtain an analytical expression for the 
initial decay rate up to second order in Xi f° r a separable 
initial state. We can also prove that the third order in % 
vanishes. Notice that for this calculation we went back 
to a general noise model with multi-body terms in G T 
(xi 7^ V I in principle). For an initial state where each 
qubit is in a pure state (/o = 1), we get 

n n 

7 = C!^( X *) 2 +c 2 J2 (4fc) 2 

3=1 k>j=l 
n 

+ c 3 £ (X-, M ) 2 + --- + °(X 4 ) (8) 

g>k>j=l 

with c v = 1 — 1/3". In JSJ we have defined the collective 
coefficients 

(x*) 2 = E (x?) 2 ; (4 fc ) 2 = E (4' fc 9 ) 2 ; etc - ( 9 ) 

and v = 1, 2, 3, . . . , n for a collective coefficient corre- 
sponding to terms with Hamming weight v. 

In the case where the coefficients are fluctuating over 
time the relevant quantity is the average 

<7)p = / l(x)PixW (10) 



, (/(* = !)) 
fo 



(7) 



since when we average the realizations of the ran- 
dom rotations, we also average realizations of the 
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fluctuating coefficients. For example for the Gaus- 
sian distribution we described before we would have 
f(Xj) 2 P(x)dx = 3(a 2 + o- 2 ), etc. Therefore for any 
distribution P the same equations hold, with the 
collective coefficients properly replaced by the strength 
parameters characterizing P. 



Eq. (JHJ) shows a decay rate 7 that is a weighted sum of 
the collective coefficients of G T . However, we are rather 
interested in obtaining a characterization distinguishing 
these coefficients. If we chose other initial states, the 
weights of the collective coefficients change; moreover, 
some vanish if some qubits are initially in the maximally 
mixed state 1/2. To make use of this feature efficiently, 
we can calculate the fidelity of the state of just a few 
qubits. Let's call M the set of m qubits that is going to 
be measured (to < n), and M its complementary (Fig. 
EJ. Thus we have 



(/< M >(*)) = Tr M [pW{t)pi M) ] 



(11) 



where we denote the reduced density matrices by = 
Trj[p]. Correspondingly, we denote as 7^ M ) the initial 
decay rate of (/< M >(t)). 

It can be shown that ^ M > is independent of the initial 
state of the qubits not being measured. This is indeed 
a desirable feature since we then don't have to worry 
about experimentally initializing them -as long as the 
separability of the initial state of the to qubits in M'^" 1 ^ 
is guaranteed. We show below the results for measuring 
the coefficients of arbitrary sets of one, two and three 
qubits, which we have labeled a, b and c; these qubits 
are initially in an arbitrary pure state. More general 
formulae is given in Appendix C; the following will suffice 
to set the basis for our proposal. 
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FIG. 5: Circuit representation of the algorithm in its final, 
experimental version. 



Po ® Po ® Po = \Va){Va \ ® Wh){Vh\ ® \<Pc){Vc\ 

i (aM = § (ixir + ixif + (x c ) 2 + 

j^a,b,c 

+ E oc-,fc) 2 + (xi jlk f + (xc tj ,k) 2 + ■ ■ 
k>j 

j.k^a,b,c 



+ g | (x:, b ) 2 + (x:, c ) 2 + (xi, c ) 2 



+ E (Xtbj) 2 + (X*a, c ,j) 2 + (Xtcj) 2 + 
j=jta,b,c 

26 



+ 27 (xlAcY + Oix 4 ) 



(14) 



It's expected that for small errors (or equivalently, for 
r small enough) only terms with low Hamming weight 
will be present in G T . This being the case, by measuring 
the initial decay rate of a few qubits, the value of se- 
lected coefficients of G T can be extracted. For example, 
if terms with Hamming weight > 3 are negligible, the 
combination 



7 « + 7 w 



7 



(a,b) 



ixU? 



(15) 



allows us to establish whether any two-body term be- 
tween an arbitrary pair of qubits a and b is present in G T . 
Notice that the measurements return the value of a given 
coefficient averaged over the distribution P{x V ji Xj'k) ( re ~ 
fer to eq. JTJJJ), giving thus its strength according to the 
parameters of P. In this way we can probe any two-qubit 
collective coefficient we arc interested in, or conduct a fair 
sampling of some of them. 

The systematic protocol to measure the collective co- 
efficients for one and two-body terms is the following: 
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1. Measure the n initial decay rates 7W) for all the 
qubits individually. That is: apply one step of 
the algorithm given in Fig. measuring only 
qubit j initially in a pure state, thus obtaining 
{f^\t = 1)). From this the initial decay rate can 
be obtained by subtracting /q = 1. 

2. Measure the n(n — l)/2 initial decay rates for all 
the possible pairs, 7(^*0. This is just as explained 
above but now measuring qubits j and k, each ini- 
tially in a pure state. 

3. With this data, all the two-body coefficients can be 
determined using i|15|) . 

4. All the one-body coefficients can be extracted by 
subtracting the two-body coefficients from the ini- 
tial decay rates of one qubit, according to |Q . 

This implementation does not distinguish between 
different product operators Oi for a given subset of the 
n qubits (i.e., between X, Y, and Z directions), since 
all the corresponding coefficients add up to form the 
collective coefficients as expressed in (0. 

If terms with higher Hamming weight are present in 
the generator, we can extend the method, but of course 
the number of initial decay rates required to map out 
the Xi increases, eventually becoming exponential in n. 
The advantage of this approach is that, when higher or- 
der terms are negligible, it makes good use of this fact. 
In addition, it provides a procedure to measure selected 
coefficients instead of going necessarily for the whole set. 
It is possible, for example, to probe the importance of 
three-body terms in G T (neglecting terms with Hamming 
weight > 4) with the combination 

7 ( a ) +7 (6) +7 (c)_ 7 ( a ,6)_ 7 (a,c)_ 7 ( 6 ,c) +7 ( ,6,c) = A^J? 

More details on the analytical calculation of the initial 
decay 7 ( M ) is given in Appendix C. For this we employed 
the tools already mentioned at the end of Sec. II. 



IV. IMPLEMENTATION ISSUES 

Consider the protocol described in Sec. Ill running on 
a canonical quantum computer: when we measure the 
state of the quantum register at the end of the computa- 
tion, we get only a binary number as a result. If the initial 
state is pure so pg M ^ = |^)(<^|, and \ip) = ® - eA/ \<Pj) with 
\<Pj) = |0) or |1), the fidelity can be evaluated in a sim- 
ple way by measuring the final state in the computational 
basis. 

The average can be implemented as follows: let X be 
the random variable which is equal to 1 if the result of the 
measurement of p^fr{t) gives Pq and otherwise. The 
random variables X u (each respresenting the outcome of 



the u th algorithm realization) all have the same distribu- 
tion Pr(X u = 1) = E(X U ) = (/ (M) ), independently of 
the particular realization. The mean 

u— 1 

(where Nr is the total number of realizations) thus gives 
an estimation of the average fidelity. From the Chernoff 
inequality we can get the minimum number of times Nr 
we must run the algorithm to achieve a precision S on 
the measure of the fidelity with an error probability less 
than e: 

Pr(|X- (/)| > S) < 2e~ 2N " s2 < e (17) 
^N R >-^ln £ - (18) 

In the scenario where three-body and higher order terms 
are neglected, we see from equation ljT5|) that to be able to 
distinguish a two-body coefficient of magnitude [3 from 
0, the number of experiments we need to run is of the 
order of lne//3 , which doesn't depend on the number of 
qubits involved and ensures the method's scalability. 

Notice again that the average has a double function. 
It's not only the average over the random rotations but 
also, when the noise is incoherent, the average of the 
fluctuating coefficients xi- 

Eqs. (|17|l and l|18|) also hold in the case where X is 
a continuous variable in the interval [ 0, 1] , in this case 
known also as the Hoeffding inequality [2(j ■ This should 
be noticed when trying to implement this protocol with 
ensemble measurements instead of projective measure- 
ments, as in NMR QIP. 

V. CONCLUSIONS 

As described in |l5j| . fault-tolerant quantum comput- 
ing requires the magnitude of the noise affecting the 
implementation of a gate to be smaller than a certain 
critical value. The quantity measuring the noise magni- 
tude and its threshold value depend on the structure of 
the noise, where by "structure" we mean which multi- 
body terms are negligible and how this scales with the 
number of qubits. Fault-tolerant thresholds of this type 
are formulated in terms of a Hamilitonian Ti. responsible 
for the errors in the computation, acting for the time 
to required to implement a gate. The generator Tit^/h 
includes the interaction with an external environment (it 
thus generates both unitary and non-unitary errors) , and 
has support on the system's space M'jss as well as outside 
of it (the environment's space). In practice, however, 
we expect to have access only to the system's space; 
the intention of our approach is actually to characterize 
a generator Qn resulting from the action of Tt in the 
system of n qubits. 
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We have presented here a protocol to analyze the noise 
structure with at most two-body terms in the generator. 
The method can be extended if higher order multi-body 
terms are present, at the price of compromising its scal- 
ability. Notice that in any case the method offers a way 
to probe the importance of these higher order terms. 

We believe our choice of the identity operator as 
primordial gate will give a fair idea of the terms present 
in the generator, keeping in mind that implemen ting 
n means implementing a time-suspension sequence [21| 
that in principle is composed of several gates modulating 
the internal Hamiltonian of the system. If we want to 
analyze the structure of the noise resulting from the 
implementation of a particular gate U g , we can easily 
account for this by implementing Ea.sE — U g U\ - its im- 
plementation will be as perfect as our ability to reverse it. 

The main advantage of our proposal is its scalability, 
but another rather important feature is that its outcome 
gives directly the generator G T of the errors over a 
time r, close to the TL referred to in the fault-tolerance 
analysis. This is one step ahead of the QPT approach, 
whose outcome is essentially a superoperator, from 
which TL still has to be extracted. At this point we 
should mention that this direct link between measurable 
quantities and the terms in an effective generator Qn 
resulting from the interaction with an environment is 
well known in the theory of relaxation in NMR systems, 
where the 1 /T2 initial decay rates encode the magnitude 
of these terms and the correlations between them 1221 ■ 



However, as we have mentioned before, our proposal 
relies on certain assumptions which limit its reach. One 
assumption is that terms with high Hamming weight can 
be neglected; this is nevertheless a reasonable one. Our 
method relies on this fact in order to achieve scalability, 
differing from QPT in that it can make good use of this 
assumption. Another assumption is that we have taken 
the error operators E to be unitary, thus confining the 
non-unitary errors we considered to a subset of unital 
processes. This assumption is again reasonable consid- 
ering the time scale of typical non-unital processes (the 
so-called Ti-processes or relaxation); unital processes are 
expected to occur faster thus the corresponding coeffi- 
cients are expected to be larger. 

Our current work is directed to generalize our proposal 
to include more general noise scenarios; also, further 
studies should concentrate in accessing the information 
we need from Ti (the full set of interaction terms for 
system plus environment) when only the system's space 
is available to us. 

This work was supported in part by the National Secu- 
rity Agency (NSA) under Army Research Office (ARO) 
contract number W911NF-05- 1-0469, by the Air Force 
Office of Scientific Research and by a Lavoisier Fellow- 
ship from the French Ministry of Foreign Affairs. 



VI. APPENDIX A: SIMPLE MODEL OF 
INCOHERENT NOISE 

Here we present an illustrative interpretation of our 
noise scenarios. The case where E is constant (coherent 
noise) evidently reflects the unitary errors that arise due 
to the implementation of imperfect gates. 

The incoherent cases we have introduced deserve some 
more analysis. A simple picture of these processes can 
be obtained by studying the effect of E for only one 
qubit, a process that is analytically tractable, being 
E = exp(— ixi<Tz)- In order to reflect the incoherent na- 
ture of the process, we of course study the evolution of the 
state averaged over different realizations of xi- Also, to 
distinguish between incoherent with long (IL) and short 
(IS) correlation times, we must observe the system at a 
time t for which the IL noise would remain constant, but 
nevertheless the IS noise will vary in each step. For these 
two processes, we have 

/oo 
(E( Xl )Y p(0) (E\xi)Y P(xi) d X i (19) 
-OO 

rOQ rOC 

Pis(t) = ■■■ E( Xl . t )...E( Xhl )p(0)EH Xhl )x 



OO J — OO 



••■ E\ Xl , t )P{ X is)...P{ X i,t) dxi,i.»dXM (20) 

where P is, as mentioned before, a Gaussian distribution 
centered in a with a deviation a a . 

If we start with a general initial state p(0) = (I + 
c x Qx + OyQy + ^ z g z )/2, we obtain 



P(t) = - + -Q Z 



-^{Qx cos(2ai) - g y sin(2a<)) 



+ -y (fe sin(2ai) + g y cos(2at)) exp(-£) (21) 

where 5 = 2a 2 t t 2 for the IL case, and S — 2a 2 t t for the IS 
one. In both cases we obtain the same physical process: 
an exponential decrease of the transversal polarization 
(together with a rotation), leaving the longitudinal po- 
larization unchanged. The decay is slower with IS errors, 
since their rapid change make them more random and 
overall less harmful. This non-unitary process has the 
following Kraus representation 



p(t) = M lP (Q)Ml + M 2 p(0)M 2 f 



Mi 



M 2 = 



l + e 



-s\ 1/2 



I e 



1 - e 



-8\ V2 



(22) 

(23) 
(24) 



This clearly shows that the process is composed of two 
parts: on one hand, a rotation (unitary operation) 
around the z axis, related to the fact that the values of 
Xi are centered around a. On the other hand, a phase- 
flip channel with probability (1 — e~ s )/2 < 0.5 of flipping 
the qubits. 
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This kind of process is commonly encountered in liquid 
NMR QIP, where this simple model reflects the essentials 
of the errors arising from spurious inhomogeneities in the 
magnetic field. In the stochastic limit, the phase-flip is 
expected to happen with a constant probability. 



VII. APPENDIX B: THE FIDELITY DECAY 
A. Analytical results 

Analytical expressions of (f(t)) can be obtained if we 
consider G T with one-body terms only, 



E t = 



E^=exp{-i X *a^) (25) 



(the actual directions hj are irrelevant) and a separable 
initial state 



Po = P^ ® Po ] ® ■ ■ • 8 P { n) 



(26) 



Under these conditions the fidelity of the whole system 
is just a multiplication of the fidelities for each qubit: 



(/(«)> =n</ w) (*)> 

We observe that 

= /W(i?p ) ,...,^ ) )d^ ) ...di?(- 



(27) 



where each integral is an average over the normalized 
Haar measure on U(2). f^\..., , . . .) is a polynomial 
function of Rs (and Rs of course). A method for 
computing this kind of averages in U (N) is presented in 
[Hill; here we will limit ourselves to state the following 
particular results 

(Tr[ARBRt]) = |Tr[A]Tr[.B] (28) 
{Tr\pR)ARpR)BR}) = \Tt[AB] (l - ^p) 

+|Tr[A]Tr[S] (Tr[p 2 ] - ±) (29) 

where all the operators belong to and we have used 
Tr[p] = 1. Applying this formulae we have 



(/»(*)) = (/»(*-!)) 



0)112 



1 



|Tr 3 [4 3) ]| 2 
6 



(30) 



This shows that f(t) only depends on the fidelity at a 
previous time t — 1, thus giving an intrinsic exponential- 
type decay. At the same time, this also shows that the 
precise decay law won't be a simple exponential. Even for 



coherent noise, where Tx^-Ej"^] is the same at all times, 
cq. <|27[) already indicates that we'll have a product of 
exponentials. 

From eq. i|5U|) it's possible to compute a closed ex- 
pression for (f(t)) in several cases, accounting also for 
the time variation of the coefficients in E t given by 
P. For coherent errors with a constant distribution 

P(Xj) = $(Xj ~ a ) we nave 



(f U) (t)) 



1 



/o 0) - § ) exp(-At) 



1 



(31) 



where 



A = - In 



4 cos 2 (a) - 1 



4a 2 



0(a 4 ) (32) 



and ffi = Tr j[p^']. To be precise, a real A like ^ is 
valid for a < tt/6 w 0.52, otherwise (f^\t)) oscillates. 
For coherent errors with non-constant distributions the 
result is very similar: for each qubit, the A in l|31|l must 
be replaced by the respective Aj . 

We also have closed expressions for the incoherent sce- 
narios proposed, including the additional averages over 
the Gaussian distribution of coefficients and already as- 
suming a,cr a << 1. For incoherent noise with long cor- 
relation time: 



(f u Ht)) 



i 



exp(— a(t)) 



1 



3°a b 



(33) 



with 



a(t) 



la 2 t 

1 + Ht 



For incoherent noise with short correlation time: 
(/ W W)«^ + (/o (i) -0exp(-^) 

with 

, /l + 2cos(2a)exp(-2cr 2 )' 
r? = — in 



(34) 



(35) 



(36) 



Fig. shows some examples of numerical calculations 
together with the theoretical result, exhibiting perfect 
agreement. 



B. Numerical results 

When multi-body terms are present in G Tl the non- 
separability of the E t prevents us from getting closed 
results for (f(t)). We studied numerically the case when 
only one-body and two-body terms are present, obtaining 
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time t (steps) 



40 60 80 

time t (steps) 



100 



FIG. 6: Exponential-like decay of (f(t)) with only one-body 

terms in G T , following our theoretical results (solid lines). For 

8 qubits initially all in the |0) state. Also we have plotted the 

numerical calculation using 100 realizations. 

□ : type IL, P Gaussian with a — 0, a a — 0.08. 

• : type IS, P Gaussian with a — 0, a a — 0.08. 

O: type C, P constant with a = 0.08. 

0: type IS, P Gaussian with a = 0.08, a a = 0.04. 

*: type IL, P Gaussian with a = 0.08, a a = 0.08. 



always a linear initial decay and a saturation value for 
long times, as mentioned in Sec. II. Figs. [7| and [S] show 
some examples of the numerical calculation of (f(t)) and 
a curve-fitting following the formula 

(/(*)> = (/o - ~) + ~ 07) 

where T is the only fitting parameter. As before, 
/o = Ti'[pq]. We don't expect an exact agreement with 
this formula; we chose it as it's the simplest exponential- 
type decay to interpolate the initial and long time behav- 
ior. Also, this choice is motivated by the fact that this 
is the exact expression for the fidelity decay when the 
random operators are rotations on U(N) (c.f. eq. (23) 
in!). 

In practice, the fitting is simply a linear fitting of 
log((/(t)) — Of course we must be careful with 

the values of (/(f)) close to (since in this range numeri- 
cally we will have null and negative values of the fidelity); 
thus we just use the points with / higher than a certain 
cut-off value f co . Notice that for fast decays, f co should 
be low enough to include a sufficiently large number of 
points to fit. 

Fig. shows the decay rate V of formula l|37|l as a 
function of the strength x, for different noise scenarios. 
It's not surprising to find that the decay is faster for 
higher x- Notice also the different proportionality be- 
tween r and x f° r the different cases. It's worth noticing 
here that the source for the data T vs. x 111 Fig- ED is 



FIG. 7: Examples of the numerical calculation of (f(t)): we 
examined the different P distributions assuming that only 
one-body and first neighbor couplings are present in G T . For 
8 qubits, initially all in the |0) state. 100 realizations. The 
fitting corresponds to formula 13711 with f co — 0.1 (see text). 

A: type C, P constant, a = 0.05, f3 = 0.02. 
□ : type C, P constant, a = 0.02, (5 = 0.05. 
O: type C, P constant, a = p = 0.05. 

*: type IS, P Gaussian with a = (3 = 0.05, a a — ap — 0.005. 
• : type IS, P Gaussian with a = f3 = 0, a a — up — 0.05. 
0: type IL, P Gaussian with a = /3 = 0.05, o a — crp — 0.005. 



the same as for 7 vs. x m Fig. El The slight difference 
lies on the nature of the fitting: in the former, we used 
eq. ()37JI with f co — 0.1 while for the later we just used 
a linear fitting of the points with / > fu m = 0.9. The 
work with T originally encouraged further studies on the 
decay rate, while later 7 became a more fruitful quantity 
for noise characterization. 



VIII. APPENDIX C: THE INITIAL DECAY 
RATE 7 

This appendix is intended to be a short guide to re- 
produce the formulas JHJ, (JEJ, fT^I . fl"4"|l and beyond, 
giving the reader more general formulae including any 
multi-body terms and also to account for the effect of a 
non-hermitian G T (complex coefficients xu cf- e( F 0)- 
We excluded this level of generality from the main text 
since the expressions are quite cumbersome. 

We work on the fidelity at t = 1, that is after only 
one iteration of the algorithm, and for a subset of m 
qubits in J*f? ( M > of the n qubits in J^j conforming the 
system under study. These are the m qubits that will be 
measured (refer to Fig. [5J). The initial state of the system 

is pg M ' ® Pq M \ where = ® jeM is separable but 
Pq M ^ could be any state. Following eqs. J3J and i|ll|) . we 
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= 0.0073 (11%) 




r 
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FIG. 8: Example of the numerical calculation of (f(t)): we 
chose a particular P distribution and we worked with different 
noise strengths. P corresponds to IL noise with coefficients 
drawn from a Gaussian distribution centered around 0. For 
8 qubits; initially, qubits 1 and 2 are in the 1/2 state while 
the rest is in the |0) state. 100 realizations. The fitting cor- 
responds to formula 13711 with f co = 0.05 (see text); the inset 
shows the obtained values of F with their relative error. 

• : one-body and first-neighbor terms in G T , a a = 073 = 0.03. 
*: one-body and first-neighbor terms in G T , a a = 073 = 0.05. 
O: with all the terms in G T , a a — 073 — 0.03. 
▲ : with all the terms in G T , a a = 073 — 0.05. 



FIG. 9: Values of F fitting the numerical data according to 
formula 13 71 . We worked with 8 qubits initially all in the |0) 
state. 

O: type C, P constant with a — j3 = x\ au the terms in G T . 
• : type IL, P Gaussian with a — /3 = 0, a a — 073 — \\ au the 
terms in G T - 

*: idem O, but only one-body and first-neighbor terms. 
□ : idem •, but only one-body and first-neighbor terms. 
<C>: idem O, but only one-body terms. 
A: idem •, but only one-body terms. 



want to calculate 

(/ (M) (* = D> = (Tr M RM) (p( E M) ) 



M. 



/W(1-7 (M) ) 



M 

(38) 



where we have made explicit the average over the random 
rotations in M and the one over the random rotations in 
M. We have used the notation <g) jeM = R^ M \ 



(M) 

Pe = Tl " 



M 




(39) 



where E = exp(— iG T ), the error operator for this first 
step, can be expanded as 



E 



l- l G T ~ l -Gl+ l -Gl + 0( X i ) (40) 
2 



G T is given by eqs. 10} and (J5J, as explained in Sec. II, 
but in principle we will allow for the xi 1° be complex 
numbers. We will rely on the fact that G T can be made 
small so high order powers of \ will be negligible; this is 
possible in theory if the time r can be made arbitrarily 
small. We thus insert l|40|) in (|39|l and keep only the 
terms up to 0(x 3 )- Then we use the separability of 



Pq M \ R( m *> and the Oi to express Ij38(l as a sum of terms 

of the form {Tr[R^ R^ R^ R^ O^}) , 

(Tr [R^p ( ( p iZtO) 0\ 3) R^p ( j) RM of 0$ ] ) , 

(Tr [R&pP R^ 0[ 3) ] ) , (Tr [R^p™ fltO') o\ j) Of 0$ ] ) 

and (TT:[R^p\pR)^0\ j) 0\P]). The average over the 
random rotations of these quantities can be computed 
using (|28[) -the last three- and (|29|l -the first two. Using 

(?) 

the fact that the 0\ are either a Pauli matrix or the 

Identity operator I, and that Trfpg ] = 1, we arrive 
to a closed expression for l|38|) in terms of the purity 



of the initial state of each qubit. 



This dependence on the purity only is not suprising; we 
can picture it the Bloch sphere: the random rotation of 
the qubits erase the information about the direction of 
the polarization vector, but not about its modulus. We 
obtain: 



\jeM jeM 



0( X 3 )(41) 
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where If the \i are rea li n °t only the second term in l|41(l van- 

{ishes but also it can be proved that the terms of cubic 

(2/3) (1 - Vj/2) if 0\ j) = a Pauli matrix order 0(x 3 ) vanish as well. 

Vj if op = I 
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